Ultracold atoms interacting with a sinusoidal mode of a high Q cavity 
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We consider the interaction of two level ultracold atoms resonant with a sinusoidal mode of the 
electromagnetic field in a high Q cavity. We found that well resolved resonances appear in the 
transmission coefficients even for actual interaction and cavity parameters. The probability of 
emission of one photon and the probability of transmission of an atom, when number or coherent 
states are initially present in the cavity, are discussed. The interplay between the increasing width 
of the resonances and multi-peak steady-state photon-statistics are also studied. Furthermore, we 
compare our results with those of a constant field mode. 

PACS numbers: 42.50.-p, 32.80.-t, 42.50.Ct, 42.50.Dv 
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I. INTRODUCTION 

Recent developments in cooling of neutral atoms |^] 
have called the attention on the interaction of ultracold 
atoms with microwave cavities [^-|7|. Very cold atoms 
(kinetic energies much smaller than the atom-field inter- 
action energy ) have a peculiar behavior when interacting 
with a cavity, they are strongly reflected unless the ra- 
tio of the size of the cavity and the atomic wave length 
associated to the interaction energy are close to certain 
resonant values. For example, it has been shown [^|||] 
that ultracold atoms of mass m may be totally reflected 
by a cavity of size L containing n photons, in a con- 
stant field mode, unless the Rabi frequency ri„ is close to 
2h{jTT / LY / (2m), j = 0, 1,2..., when only approximately 
half of the atoms are reflected. 

Models for the interaction of atoms with a microwave 
cavity usually do not take into account the spatial de- 
pendence of the interaction coming from the electric field 
profile 1^. This is justified because, in general, we are in- 
terested in the atom-field state after the atom has crossed 
the cavity and are dealing with atomic thermal velocities, 
which are exceedingly high compared with the recoil ones 
experienced by the atom inside the cavity. 

In order that we may observe any effect on the center- 
of-mass motion of the atoms, their velocities have to be 
reduced to very low values in such a way that even the 
microscopic inversion of momentum of the center-of-mass 
could take place. Actually, the only way for an atom 
to experience mechanical action on its center-of-mass is 
through the mechanism of exchange of momentum with 
the whole cavity. At high atomic velocities the energy 
exchange between the field and the atoms takes place 
through the Rabi flopping of the internal atomic levels. 
At very low atomic velocities the spatial variation of the 
electromagnetic field acting on the atoms provides an ad- 
ditional driving mechanism which lies in the momentum 
exchange of the atom with the cavity. 

The aim of the present work is to investigate the effect 



of a sinusoidal electric field mode when ultracold atoms 
are sent into a cylindrical microwave cavity for actual 
experimental parameters. We will show that under these 
conditions the behavior of the relevant physical quantities 
may be quite different from those of a discontinuous field 
mode, although still preserving resonant features if the 
atoms are cold enough. 



II. THE MODEL 

Let us assume that a two level atom is travelling along 
the z-direction in the way of a microwave high Q cavity of 
length L, containing a quantized field mode in resonance 
with the two level atomic transition |e) <^ \f). 

The Hamiltonian describing the motion of the two level 
atom in the z -direction is given, in the interaction rep- 
resentation, by 
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Here a = |/)(e|, a (a^) is the annihilation (creation) op- 
erator of the electromagnetic fleld mode whose spatial 
modulation pattern is u{z) and Hk'^ / {2m) is half the Rabi 
frequency associated with the transition |/) |e) in the 
vacuum. u{z) is normalized in such a way that the area 
under the mode is L. In refs. |^,|| u{z) is taken as a 
square function of height 1 and width L (mesa function) . 
In this case it can be shown that the transmission of ul- 
tracold atoms through a cavity containing n photons is 
very small unless 



Kn = K-^n + 1 



(2) 



is a multiple of n/L. In ref. |^ the sech '^{z/L) and sinu- 
soidal modes were considered and the authors have con- 
cluded that, for high values of L, the resonances referred 
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above are smeared out. Here we will discuss in more de- 
tail the sinusoidal mode having just one antinode. We 
will show that in this case, resonances still exist for high 
values of k„L as long as the atoms are cold enough. 

The Hamiltonian given in Eq. ^ can be conveniently 
written in the dressed state basis as: 



where V±„(z) is given by 



and 



|±,n) = (|e,n)±|/,n+l))/V2. 



(3) 
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In the dressed state basis our problem is equivalent 
to that of a particle being scattered by the potentials 
V±n{z). The cavity acts as a barrier when the atom- field 
state is \ +n) and acts as a well when the atom- field state 
is \ — n). 

Before the atom reaches the cavity containing the field 
^ c{n)\n) , the initial atomic state is given by the product 
of a wave packet in momentum space (z|0)= J dk(f){k)e'^^^ 
times the internal energy state of the atom that, for 
definitiveness, we assume to be the upper state |e). After 
the atom interacts with the cavity the total atom-cavity 
state is (for \z\ > L/2): 



Vc(n) / dfc</)(fc)exp(-»:^r)[(e*^-t+„(A;) + 
\+n) + {e'^^H^^^ik) + e-''=V_„(fc))| - n)], 



Ak)) 
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where the reflection, r±n, and transmission, t±n, coeffi- 
cients are associated with the scattering eigenfunctions 
if^{z). These functions are solutions of the time inde- 
pendent Schrodinger equation : 



dz^ 



(7) 



where hk is the momentum of the incident particle. 

In this paper we are interested in discussing a situation 
where the longitudinal z-dependence of the electric field 
is sinusoidal: 



u{z) = 9{L/2 - \z\){n/2) cos( — ). 

L 



(8) 



This is the case of a cylindrical cavity tuned cither to a 
TE,„„i or a TM rnni mode ^]. 

We may try to obtain the values of r±n{k) and t±n{k 
by numerically integrating a complex solution of Eq 
for runaway asymptotic conditions. An alternative way, 
where we found numerical and analytical advantages, is 
to consider instead the even and odd (v?^"(z)) 

real (up to an overall complex constant) eigensolutions 
of Eq. 0, such that for z < -L/2 @ 



^±„"(z)«cos(fcz-5±„"). 



(9) 



Here S^^{k) ((5^"(fc)) is the phase shift associated to the 
even (odd) eigensolution of the Schrodinger equation. In 
terms of these phase shifts the reflection and transmission 
coefficients may be written as 



r±„(fc) 



t±n{k) = 



(10) 



The phase shifts are related to the logarithmic deriva- 
tives at z = —L/2 by 



f3^: = k tan(«„i/2+5±„") 



(11) 



Notice that f3^" are real in contrast with the complex 
logarithmic derivative of a typical, right or left, propa- 
gating solution. /From Eqs. 10 and ^it is easy to show 
that: 
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In the case of the potential barrier, for K„i ^ 1 and 
k <^ Km the transmission coefficient t+„ is very small 
and close to zero. In the case of the potential well, and 
for k <C Kn the transmission coefficient i„„ is close to 
zero unless /3^"/3(7" = —k"^, that is for very small /3~"or 
very small For a fixed value of k/Kn, these res- 

onances in the transmission coefficient occur at special 
values of k^L, which would correspond to the appear- 
ance of a bound state with fc^ = ( notice that the 
transmission coefficient has a pole at fc = iPeo)- 

Assume that initially there are n photons inside the 
cavity and the atom is incident from the left in the upper 
state |e). After the atom interacts with the cavity, it 
can be transmitted in the upper state, transmitted in 
the lower state, reflected in the upper state or reflected 
in the lower state. The probabilities that these events 
occur are denoted by T", TJ?, i?"and RJ, respectively 
and given by: 

T:^]\t+^+t_„\\ T;-i|t+„-t_„|2 (14) 



and 



R:^j\r+n + r^nr, R] = -\r+^ - r^n\. (15) 

Of course TJ' + + R'^ + R] = 1. We will show in 
the next section that for /c/k„ ^ 1, t+,i = 0. In this 
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case T" « Tj « |t_„p/4. Due to this last approximate 
relation and the unitarity condition we choose to discuss 
below only two quantities of physical interest: the prob- 
ability of emission of a photon, P^^^ — Tj + R^, and the 
probability of transmission of an atom, T" — + TJ . 
We notice also that when r_„ = 0, PJJ,, ^ T". 

If we initially have a photon number distribution given 
by |cnP, the total probability of photon emission is 
^cm = X] knP^'c'm ^"^^ the total probability of atomic 
transmission is T = ^ |c„pT". These quantities should 
be easily measurable by detecting the internal energy lev- 
els of the reflected and transmitted atoms. As we will 
show, these quantities depend strongly on the shape of 
the field inside the cavity, as long as fc/K„ ^1. 



III. ANALYTICAL RESULTS 



The constant field mode {u{z) 
have been discussed in refs. (2lp,H 
have « -K„ tanh(K„i/2), (3f' 
/3r" ~ tan(K„L/2) and 
that the resonances appear when 



= I, for \z\ < L/2) 
\. For k/Kn ^ 1 we 
' « -K„ coth(K„L/2), 
~ K„ cot(K„i/2), so 
KnL equals to an in- 



teger multiple of tt yj^,|5|. In this case the resonances 
are very sharp, their width being approximately constant 
and equal to 4fc/K„ even for large k„L. 

In the case of the sinusoidal mode, as long as k„L is 
small, it is easy to obtain the even and odd eigensolu- 
tions of Eq. ^ by numerical integration . This can be 
done easily and with confidence for IOOtt. If we 

wish to predict values closer to real experimental situa- 
tions we need to consider higher values of k„L, which are 
typically of the order of 0(10^ — 10^) for Rydberg atoms. 
In these cases the numerical solutions do not converge 
rapidly for the | — n) channel and we found easier, and 
more instructive, to use a WKB-like approximation for 
calculating the S^" and Sf" phase shifts or, equivalently, 
the logarithmic derivatives 

In the case of the potential barrier, and for fc/K„ ^1, 
we may use the WKB approximation |jll|] for calculating 
/3+^ . We get 



+n 



-fc(iTe), 



(16) 



where the — and -I- signs inside the parenthesis refer to 
the even and odd solution and 



Q = exp(- 



(17) 



where ±a are the turning points. Substituting the val- 
ues for in Eqs. ^ and 13 we get r+„(fc)~— ie^*'"'^ 
and t+n{k)~ iQe'^'^^ . For fc/K„ < 1, o ~ L/2 and 
Q ~exp(— 4.19k„L), which is very small for k„L > 1. 

Inside the cavity the solutions if~'2{z) for the potential 
well may be written, whenever the semiclassical approx- 
imation holds, as 



-^=cos{ q{z')dz') 



r-r^ sin( / q{z')dz'), 
\/q{z) Jo 



(18) 



where q{z) = y^k^ + k\u{z). It is well known that Eqs. 
[T^ are valid when the condition 



dq{z) 



dz 



(19) 



is fulfilled ||Tl|. This condition may be satisfied for the 
sinusoidal mode and for any z, whenever 



7r^(Kn/fc)' 



< 1. 



(20) 



In this case the semiclassical solutions (Eqs. ) are 
valid for all z and the reflection coefficient r^n{k) van- 
ishes as ^ goes to zero, as we will show below. When ^ > 1 
the inequality [l9| is not valid near the regions \z\ < L/2. 
In this case the modulus of r_„(fc) may increase even to 
1. Notice that we may have ^ > 1 even for large kL and 
KnL. For example, for k„L ~ 10^ and kL ~ 10^, ^ « 24. 

Near z = —L/2, u{z) may be approximated by a 
straight line and the solutions of Eq. [7| may be written 
as: 

(^-^(z) CX wl/3(A,o Jl/3(l«) + Se,oJ-l/3(w)) (21) 

where 

4fc2 X 3/2 



TT , / 2Z 
W = —KnL 
6 



1 



(22) 



The logarithmic derivatives at z = —L/2 are given by: 

Ae,oJ[/,{r'/i)+Be,oJL^/,{r'/-i)\ 



^e,oJl/3(^-V3) + Be,oJ-l/3(e-V3) 



(23) 



By connecting the even and odd solutions of Eq. gT| with 
those given by Eqs. |l^ , we obtain : 



k ^ ^^^^ 



sin((/7 + 7r/12) 



cos(iy9 — 7r/12) 
and Ao/Bo = xif + '''/2), where 

-L/2 



+ K'^^u{z)dz 



(24) 



(25) 



For ^ ^ 1, we use the asymptotic expressions for 
J±i/?,{C^ I'i) to obtain /S"" w -k{i/2 + ism^p) and 
« -fc(C/2 - cot Lp). From Eqs. H and [l| we obtain 
the same result (up to a phase coming from our definition 
of r_„ and i_„) of ref. 
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-i(C/4) [cos 2lp + (^/4) sin 2if] t. 



(26) 



t- 



-ikL 



[i^/4.fe~^"^ + (1 + ie/4)^e2»'^] 



For ^ 3> 1, we may use the series expansion for 
•^±1/3(^73) to obtain /J"" = -ak^^^^xi^) and = 
afc^i/3x(v3 + 71-/2) where a = (2/9)i/3r(2/3)/r(4/3). In 
this case we get: 



(1 + tae^^'xivW - + V2)) ' 



ikL 



(27) 



and 



l-a^g^/My>)x(y> + 7r/2) 



^—ikL 



(28) 



Resonances occur when or are zero, that is 
when either xiv) or xif + '^/'^) ^'^^ zero, which give us 
a simple condition for locahzing them: 



K„,L 



7r/2 



^/{k/Kny + {Tr/2)cos{d)d9 = m7r/2 + 7r/12, 

(29) 



with m = 0,1,2,... For {k/K„)^ < 1 we may evaluate 
the integral approximately and get for the resonance po- 
sitions: 



KnL - 2.092 X (TO7r/2 + 7r/12). 



(30) 



Therefore, for the sinusoidal field mode we will also 
obtain a series of resonances at values of that are 
separated by an interval approximately equal to tt. In 
addition their positions are shifted in relation to the con- 
stant field case, their widths are wider and increase with 
KnL. We will show that these results give rise to new 
features in the physical quantities of interest. 



IV. NUMERICAL RESULTS 

We have checked that the WKB-like eigenfunctions 
agrees extremely well with the numerical solutions at 
each point, even for low values of k„L and k/Kn- For 
higher values of k„L (realistic values included) the WKB- 
like solution should be increasingly better. 

For very small values of k„L the resonances in , 
although wider than in the constant field case, do not 
correspond to appreciable values of for m ^ n. For 
large values of k„L the resonances in are wide enough 
to produce appreciable values of for m ^ n , so one 
should take them in account when predicting the total 
probabilities Pcm , or T . This is not the case for the 



constant field mode where the narrow resonances assure 
a similar behavior of PJ^ and Pem, and T, except for 
accidental numerical coincidences. 

In Fig. 1 we show our results for the probability of 
emission of a photon when there is zero photons inside 
the cavity (vacuum field), P°^^ for a fixed value oik/ Kn — 
0.01, as a function of k„L in the range IOOtt < k„L < 
1047r. We see that the resonances are shifted and wider 
in the sinusoidal mode (solid line), in comparison with 
the constant mode (dashed line), and that their widths 
are still much smaller than the resonance separation. 

We now consider that the cavity is in contact with a 
heat bath at a temperature associated with a thermal 
mean photon number rtf,. Assuming, as usual, that we 
may treat gain and loss independently and that the inter- 
val of the atom separation obeys a Poissonian statistics, 
it can be shown that the stationary photon distribution 
is given by [| 



dpn 
dt 



r(v iP""-^ - V P" 

' V-Fn— 1^ cm Pn^ cm 



:K + 1) 



[npn - {n + l)p„+i] - Q'Tibiin + l)p„ - npn+i]. (31) 

Here r is the atomic injection rate, Q is the cavity quality 
factor, Lu is the cavity resonance frequency and n?, is the 
mean number of photons inside a cavity in thermal equi- 
librium with a reservoir at temperature T^. From Eq. ^ 
we obtain the stationary population : 



Pn = Po 



n 



{nb + 1) 



(32) 



where Ni^^ — rQ/uj. This result is similar to that ob- 
tained in the conventional micromaser, where the value 
of PJJ^ is given by sin^(rK^/(2m)), r — {m/hk)L being 
the time of flight of the atom through the cavity. For 
fc/Kp <^ 1 the photon distribution obtained from Eq. 
p3 is completely different from that of the conventional 
micromaser. In this case the predictions are strongly de- 
pendent on the field mode profile as long as Ub is small 
enough and N^x large enough. 

Fig. 2 shows the stationary photon distributions, for 
Nex = 1000, rif, = 1 and /c/kq = 0.01, when the K2L 
value corresponds to the position of the 100*'' resonance 
in the case of the constant mode {K2L = 99tt) and of 
the sinusoidal mode {K2L sa 103. Ttt). Fig 2a (constant 
mode) shows the apparition of two thermal distributions 
with peaks at n — 3,12. For such low values of kL 
this may happen in two special cases when there 
is a resonance as in nL^fi = (997r/3^/^)(3)^/^ = 997r 
{n — 3) or when there is a numerical coincidence as in 
kL^/VI = (99^31/-*) (12)1/4 ^ i40.00 77r (n = 12). Oth- 
erwise we would obtain an equilibrium thermal distri- 
bution at the temperature corresponding to the average 
photon number n^. This is not the case for the sinusoidal 
mode (see Fig. 2b) since the probability of photon emis- 
sion, P^j^ , may be important for several different values 
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of j besides the resonant probability for n = 0, even at 
these low values of kL. This is due to the large width of 
the resonances for the sinusoidal mode, in contrast with 
the sharp ones for the constant field case. As we men- 
tioned before, these widths increase with kL for the si- 
nusoidal mode, producing an extremely complex photon 
statistics. 

Realistic values for the coupling of Rydberg atoms with 
microwave cavities are much larger than the values con- 
sidered above. For example, consider the Rydberg transi- 
tion between circular states with principal quantum num- 
bers 50 and 51 in Rubidium at 51GHz, used by the Ecole 
Normale Superieure Group in recent experiments [ p2| . 
For a cylindrical cavity, with a length L w QAOcm and a 
radius R w 0.68cto, this transition is resonant with the 
TEi2imode. For a dipole moment d w IQ^^^Cm, we get 
Ko ~ 2.1 X lO^cm^iand kqL w 270007r. 

In Fig. 3 we plot the probability of detecting one trans- 
mitted atom, T" = + Tf, and the probability that 
one photon be emitted, Pcm= + TJ, if we have ini- 
tially n = 0,1,2,3 in the cavity, as a function of KqL, 
for 300007r < kqL < 300057r and k/no = 0.01. This ra- 
tio corresponds to Rubidium atom velocities of approxi- 
mately 0.2mm/s. Atomic velocities lower than that have 
been already obtained by evaporative cooling Q. For 
these velocity values the injection rate should be very 
small, if we wish to have only one atom at a time in 
the cavity (r ~ 0.06s^^), and the value of Nex would 
be smaller than 1 even for Q « 3x10^^! In this case the 
stationary photon distribution would be very close to a 
thermal one, inhibiting the new features shown in Fig. 
2. In Fig. 3, after a careful inspection, we note that 
T" and show a similar, but not identical, behavior. 

presents better resolved resonances than T" and this 
could be useful for determining the more convenient ex- 
perimental setup. 

In Fig. 4 we plot the total probability of detecting 
the transmitted atom, T, for the sinusoidal mode (dotted 
line) and the constant mode (solid line), when a coherent 
state, either with n = 0.25 or n = 2, is initially present 
inside the cavity. The figure shows very sharp resonances 
in KqL for the constant mode. They are a consequence 
of the fact that only one resonance is contributing for a 
given n, with a statistical weight |c(n)p. For the sinu- 
soidal mode the resonances in T" are broader and several 
resonances may contribute for a fixed kqL. The case in 
which n = 2 is much more sensible to these features, 
showing, for the constant mode, many little resonance 
peaks and, for the sinusoidal mode, very bad resolved 
resonances. 



For warm atoms it is the value of the Rabi angle and the 
rate of incident atoms that defines the gain. For ultracold 
atoms we should also consider the exchange of momen- 
tum between the cavity and the atom. Most atoms are 
reflected when the atom- field system is in the state \+n), 
independently of the shape of the mode. When the atom- 
field system is in the state | — n) and the shape of the 
cavity is smooth enough, so that the effective de Broglie 
wavelength does not vary appreciably through the cavity, 
most atoms would be transmitted. When this is not the 
case the atoms would be reflected, unless close to 

certain resonant values. 

We have calculated, in the case of a sinusoidal mode, 
the probability of transmission of an ultracold atom and 
the probability of emission of a photon comparing them 
with the corresponding mesa function results. In the case 
of a mesa function sharp resonances appear in the trans- 
mission coefficient as a function of k„L. In the case of the 
sinusoidal mode resonances are present when the value of 

4^K^L^ > 1, even when k„L is very large. This may 
give us a chance to test these effects in a realistic ex- 
periment (when KqL ~ 10^) using a cylindrical cavity. 
We also found, in the case of the sinusoidal mode, that 
the width of the resonances are larger, but still well re- 
solved for realistic parameters, and increases with kqL, 
as long as k/Ko is very small. As a consequence, the 
transmission of atoms through a cavity containing a co- 
herent state with few photons, as a function of kqL, has 
a completely different shape for the two modes we have 
studied. The broadening of the resonances produces also 
a more complex photon statistics for the stationary state, 
at least for theoretical parameters, to our knowledge, far 
from realistic ones. 

These effects show the importance of considering the 
details of the cavity mode when analyzing the scattering 
of ultracold atoms by a resonant cavity. 

We should remark that a serious difficulty in realizing 
these experiments is the long time of flight that ultracold 
atoms need to reach the cavity after their production. 

VI. ACKNOWLEDGMENTS 

This work is supported in part by Centro Latino Amer- 
icano de Fi'sica (CLAF), The Brazilian Conselho Nacional 
de Desenvolvimento Cientffico e Tecnologico (CNPq), 
Programa de Apoio a Niicleos de Excelencia (PRO NEX) 
and Fundagao Universitaria Jose Bonifacio (FUJB). 



V. CONCLUSIONS 



The behavior of ultracold atoms interacting with a field 
inside a high quality microwave cavity depends in an es- 
sential way on the profile of the mode inside the cavity. 
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FIG. 1. Probability of emission of a photon when the cavity 
is initially empty as a function of kL when the mode is: the 
sinusoidal mode (solid line) and the constant mode (dotted 
line) 

FIG. 2. Stationary photon distribution for Nex = 1000, 
k/K = 0.01 and = 1 when K2L corresponds to the 100*'' 
resonance: (a)constant mode; (b) sinusoidal mode. 

FIG. 3. (a) Probability that an atom being transmitted 
through the cavity and (b)Probability of emission of a photon 
when the cavity has initially n photons, as a function of kL: 
n = (solid line), n = 1 (dashed line), n = 2 (dotted line). 
k/K = 0.01. Sinusoidal mode 

FIG. 4. Probability that an atom being transmitted when 
initially there is a coherent state in the cavity as a func- 
tion of kL for the constant mode (full line) and the cosine 
mode (dashed line ) (a) : n = 0.25 (solid line); (b): n = 2. 
k/K = 0.01. Sinusoidal mode 
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